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ABSTRACT 

We report the discovery of three new transiting extrasolar planets orbiting moderately bright (V = 
11.1 to 12.4) F stars. The planets have periods of P = 2.6940d to 4.4572d, masses of 0.60Mj 
to 0.80Mj, and radii of 1.57i?j to 1.73 They orbit stars with masses between 1.40 Mq and 
1.51 Mq. The three planets are members of an emerging population of highly inflated Jupiters with 
0.4 Mj < M < 1.5 Mj and R > 1.5 

Subject headings: planetary systems — stars: individual (HAT-P-39, GSC 1364-01424, HAT-P-40, 
GSC 3607-01028, HAT-P-41, GSC 0488-02442) — techniques: spectroscopic, pho- 
tometric 



1. INTRODUCTION 

Transiting exoplanets (TEPs) are key objects for the 
study of planets outside the Solar System. The geom- 
etry of these planetary systems enables measurements 
of several important physical parameters, such as plane- 
tary masses and radii, or the sky-projected angle between 
the orbit al axis of a planet and the spin axis of its host 
star (e.g. iQueloz et al.ll2006T ). The vast majority of well- 
characterized TEPs (i.e. TEPs with measured masses 
and radii) have been discovered by dedicated photomet- 
ric surveys, including the Wide Angle Search for Planets 
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(WASP: iPoilacco et al.|[2006D . the Hu ngarian -made Au- 
tomated Telescope Network (HATNet: iBakos et al.ll2004D 
and its sout hern extension (H ATSouth; Bako s et al 
20l2]) . Kepler dBorucki et aljf2010h CoRoT (Barge et al 



2008)^ OGL E (lUdalski et al.ll2002h. T rES (jAlonso et al 
2004), XO (iMcCullough et alJ I2006D. the Qatar Exo 



planet Survey (QES: IAlsubai et al.ll201lD . t he Kilodegree 
Extre mely Little Te lescope survey (KELT: [Siverd et al.l 
HoHI), and MEarth (jCharbonneau et al.ll2009ir ~ 

Significant among these are the ground-based, wide- 
field surveys using small aperture telescopes, includ- 
ing WASP, HATNet, HATSouth, TrES, XO, QES, and 
KELT. While these surveys are heavily biased towards 
discovering large planets on short-period orbits com- 
pared to the Kepler and CoRoT space-based surveys, the 
planets discovered by ground-based surveys tend to or- 
bit stars that are brighter than those discovered by the 
space-based surveys, making such planets more amenable 
to detailed characterization and follow-up studies (this 
is true as well for the few, but valuable, transiting plan- 
ets discovered by RV searches, which are found around 
even brighter stars than those discovered by photomet- 
ric surveys). Additionally the extreme environments in 
which these planets are discovered, while perhaps not 
representative of most planetary systems, create a nat- 
ural experiment for testing theories of planet structure 
and formation. For example, a number of gas-giant plan- 
ets have been discovered with radii tha t are substantially 
large r than theoretically expected (e . g.lMandushev et al.l 



20071: iCollier Cameron et al.l 120071: iSnellen et al.l 120091: 



Hebb et al.l 120091: lLatham et aU 120101: iFortney et ail 



201U lAnderson et all I2010L l201lMEnoc 
Hartman et al.l 120111: ISmallev et all 120121) 



Enoch et all 120111: 

These have 

been used to empiricall y determine the fa ctors affecting 
the radii of planets (e.g. lEnoch et al.ll2012f ) which in turn 
informs theoretical work on the subject. 

In this paper we present the discovery and character- 
ization of three new transiting planets around the rel- 
atively bright stars GSC 1364-01424, GSC 3607-01028, 
and GSC 0488-02442, by the HATNet survey. As mem- 
bers of the growing sample of highly inflated planets, 
these objects will provide valuable leverage for under- 
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standing the physics that determines the structure of 
planets. 

In Section [2] we summarize the detection of the pho- 
tometric transit signal and the subsequent spectroscopic 
and photometric observations of each star to confirm the 
planets. In Section |3] we analyze the data to rule out 
false positive scenarios, and to determine the stellar and 
planetary parameters. Our findings arc briefly discussed 
in Section |4j 

2. OBSERVATIONS 

The observational procedure employed by HATNet 
to discover TEPs has been describe d in detail in sev- 
eral previous discov ery papers (e.g. iBakos et al.l 120101 : 
lLatham et al.l 120091 ) . In the following subsections we 
highlight specific details of the procedure that are rel- 
evant to the discoveries of HAT-P-39b through HAT-P- 
41b. 

2.1. Photometric detection 

Table [T] summarizes the photometric observations of 
each new planetary system, including the discovery ob- 
servations made with the HATNet system. The HATNet 
images were processed and reduced to trend- filtered light 
curves following the procedure described by IBakos et al.l 
(2010). The light curves were searched for periodic box- 
shaped signals usi ng the Box Least-Squares (BLS; see 
iKovacs et al.ll2002l ) method. We detected significant sig- 
nals in the light curves of the stars summarized below 
(see Figure [T]) : 

• HAT-P-39 - GSC 1364-01424 (also known as 
2MASS 07350197+1749482; a = 07 h 35 m 01.97s, 

5 = +17°49'48 .3"; J2000; V = 12.422, 
iDroege et al.ll2006T l. A signal was detected for this 
star with an apparent depth of ^10.9 mmag, and a 
period of P — 3.5439 days. 

• HAT-P-40 - GSC 3607-01028 (also known as 
2MASS 22220308+4527265; a = 22 h 22 m 03.00s, 

6 = +45°27'26 .6"; J2000; V = 11.699, 
IDroege et al.ll2006T ). A signal was detected for this 
star with an apparent depth of ~4.4mmag, and a 
period of P = 4.4572 days. 

• HAT-P-41 - GSC 0488-02442 (also known as 
2MASS 19491743+0440207; a = 19 h 49 m 17.40s, 
5 = +Q4°40 / 20 .7 // ; J2000; V = 11.087, 
IDroege et al.H2006T ). A signal was detected for this 
star with an apparent depth of ~8.4mmag, and a 
period of P = 2.6940 days. 

2.2. Reconnaissance Spectroscopy 

High- resolution, low-S/N "reconnaissance" spectra 
were obtained for HAT-P-39, HAT-P-40, and HAT-P-41 
using the Harvard-Smithsonian Center for Astr ophysics 
(CfA) Digital Speedometer (PS; lLatham] 119921 ) until it 
was retired in 2009, and thereafter the Tillingha st Re- 
flector Echelle Spectrograph (TRES: lFureszll2008f ) . both 
on the 1.5 m Tillinghast Reflector at the Fred Lawrence 
Whipple Observatory (FLWO) in AZ. The reconnais- 
sance spectroscopic observations and results for each sys- 
tem are summarized in Table [2] The DS observations 
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Fig. 1.— HATNet light curves of HAT-P-39 (top), HAT-P-40 
(middle) and HAT-P-41 (bottom) phase folded with the transit pe- 
riod. In each case we show two panels: the top shows the unbinned 
light curve, while the bottom shows the region zoomed-in on the 
transit, with dark filled circles for the light curve binned in phase 
with a binsize of 0.002. The solid line shows the model fit to the 
light curve. 

were redu ced and analyzed fo llowing the procedure de- 
scribed bv lTorres et~aTT (|2002[ ) . while the TRES observa- 
tions were red uced and ana l yzed followin g the procedure 
descri bed by lOuinn et al.l (|2012t) and iBuchhave et al.l 
(|2010t) . 

Based on the observations summarized in Table [2] 
we find that all 3 systems have RV root mean square 
(rms) residuals consistent with no detectable RV varia- 
tion within the precision of the measurements. All spec- 
tra were single-lined, i.e., there is no evidence that any 
of these targets consist of more than one star. Note 
that while there is a close companion to HAT-P-41 
(Section |2~5|) . it was resolved by the TRES guider and 
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TABLE 1 

Summary of photometric observations 



Instrument /Field 
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the light from the companion did not go down the fiber. 
The gravities for all of the stars indicate that none of 
the stars are giants, though HAT-P-40 may be slightly- 
evolved. 

2.3. High resolution, high S/N spectroscopy 

We proceeded with the follow-up of each candidate by 
obtaining high-resolution, high-S/N spectra to charac- 
terize the RV variations, and to refine the determination 
of the stellar parameters. Th e observations were made 
with HIRES (jVoet et al.lll994H on the Keck-I telescope in 
Hawaii, and with FIES on the Nor dic Optical Telescope 
on th e island of La Palma, Spain ; Diunvik fe; Andersen! 
120101 ). We used the high- resolution fiber (providing spec- 
tra with a resolution R — 67,000) for four of the FIES ob- 
servations, and the medium-resolution fiber (R = 46,000) 
for five of the FIES observations. The HIRES observa- 
tions were reduced to radial velocities in th e barycentric 
frame following the procedure described by IButler et al.l 
(1996), while the FIES obser vations were reduced fol- 
lowing iBuchhave et al.l (|2010f ). The RV measurements 
and uncertainties are given in Tables [5][5] for HAT-P- 
39 through HAT-P-41, respectively. The period- folded 
data, along with our best fit described below in Section [3] 
are displayed in Figures 

In each figure we show also the spectral line bisector 
spans ( BSs) comput e d from the Keck/HIRES spectra fol- 
lowing (^rreseTil] (|2007l ). the full width at half maxi- 
mum (FWHM) of the Keck/HIRES spectral lines (com- 
puted from the cross-correlation function in a similar 
manner t o the BSs), and the S acti vity index calculated 
following llsaacson fc Fischer! ()2010D . 



2.4. Photometric follow-up observations 

We conducted additional photometric observations of 
the three stars with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope, the Spectral CCD on the 2.0 m 
Faulkes Telescope North (FTN) at Halcakala Observa- 
tory in Hawaii, and the CCD imager on the Byrne Ob- 
servatory at Sedgwick (BOS) 0.8 m telescope, at Sedg- 
wick Reserve in the Santa Ynez Valley, CA. Both FTN 
and BOS are operated by the Las Cumbres Observa- 
tory Global Telescope (LCOGlB Brown et al. 2012, in 
preparation). The observations for each target are sum- 
marized in Table [1] 

The reduction of the KeplerC am images to light curves 
was performed as described by .Eakoset al. (2010). The 
FTN and BOS images were reduced in a similar manner. 
We performed EPD and TFA to remove trends simulta- 
neo usly with the light curve modeling (for more details, 
see iBakos et al . (2010)). The final time series, together 
with our best-fit transit light curve model, are shown 
in the top portion of Figures [M3 while the individual 
measurements are reported in Tables [~ 



2.5. Adaptive Optics Imaging 

We obtained high-resolution imaging of HAT-P-41 on 
the night of 21 June 2011 using the Clio2 near-IR imager 
on the MMT 6.5 m telescope in AZ. Observations were 
obtained with the adaptive optics (AO) system in H- 
band and in L'-band. Figure [S] shows the resulting H- 
band image of HAT-P-41 which easily resolves the 3'.'56± 
0'.'02 neighbor. 

16 http://lcogt.net 
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TABLE 2 

Summary of reconnaissance spectroscopy observations 



Instrument 


HJD - 2400000 


-*■ err * 




v sin % 


/rtv 






(K) 


(cgs) 


(kms- 1 ) 


(km s _1 ) 


HAT-P-39 












DS 


4 obs 54807-54931 


6250 


4.0 


16 


28.54 ± 0.58 (rms) 


TRES 


54934.6546 


6500 


1.0 


16 


29.25 


HAT-P-40 












TRES 


55084.8821 


6110 ± 80 


4.21 ± 0.13 


8.7 ± 0.5 


-25.97 


TRES 


55131.6811 


5940 ±170 


4.04 ± 0.26 


10.8 ± 1.5 


-25.67 


TRES 


55138.6609 


5910 ± 170 


3.81 ±0.27 


12.4 ± 1.6 


-26.08 


TRES 


55162.6713 


6020 ± 50 


3.93 ±0.10 


8.0 ± 0.5 


-25.65 


TRES 


55168.5775 


6120 ± 80 


4.15 ±0.14 


8.5 ±0.5 


-25.58 


HAT-P-41 












TRES 


55319.9727 


6504 ± 100 


4.3 ±0.16 


23.9 ±0.7 


32.32 


TRES 


55372.9209 


5807 ± 223 


3.94 ±0.35 


32.1 ± 2.5 


29.84 


TRES 


55373.9046 


6430 ± 105 


4.28 ±0.17 


27.5 ±0.7 


30.44 



a The mean heliocentric RV of the target on the IAU system, with a systematic uncertainty of 
approximately 0.1 kms -1 mostly limited by how well the velocities of the standard stars have been 
established. We give the mean and rms RV for the four DS observations of HAT-P-39, the velocity 
and classification for each TRES observation of HAT-P-39 through HAT-P-41 is listed individually. 



TABLE 3 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-39. 



BJD a RV b <trv c BS o-bs S d Phase 



(2,454,000+) 


(m s 


') 


(m s 


- 1 ) 


(m s 


- 1 ) 


(m s 


- 1 ) 










954.81362 


-51 


.56 


7. 


.81 










0. 


180 





.345 


955.84224 


32 


.48 


8. 


.82 










0. 


171 


0. 


.635 


956.82792 


61. 


.50 


9. 


51 










0. 


175 


0. 


.913 


1107.13759 


-90. 


.39 


8. 


88 


-4. 


.63 


4. 


91 


0. 


197 





.327 


1109.07573 


78 


.38 


9. 


.05 


-13 


.59 


7. 


.57 


0. 


183 





.874 


1112.12944 


42 


.24 


7. 


.88 


-4 


.74 


5. 


.19 


0. 


176 





.736 


1135.07184 


-39 


.87 


6. 


.92 


-2 


.61 


4. 


.96 


0. 


178 





.210 


1172.95852 


-41 


.57 


7. 


.30 


19. 


77 


1. 


.58 


0. 


176 
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1187.99450 


-28. 


.34 


8. 


.69 


8. 


.74 


5. 


.51 









.143 


1190.10225 


70 


.06 


7. 


12 


30. 


.09 


4. 


.51 


0. 


181 





.738 


1192.03996 


-146 


.19 


6. 


26 


16. 


.58 


6. 


.61 


0. 


181 





.285 


1192.05453 










16. 


.75 


3. 


.03 


0. 


180 





.289 


1193.05376 


-31 


.04 


6. 


.96 


-r 


.15 


4 


07 


0. 


179 





.571 


1193.88084 


90 


.74 


7. 


.48 


-36 


.37 


5. 


.72 


0. 


178 





.804 


1197.00264 


89 


.62 


7. 


.37 










0. 


180 





.685 


1198.91700 


4. 


.31 


10. 


.58 










0. 


184 





.225 


1250.89571 


129. 


15 


7. 


.46 


-8. 


.86 


4 


.93 


0. 


183 





.892 


1251.90337 


-98 


.98 


7. 


.98 


-6 


.43 


4 


.96 


0. 


179 
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1289.75983 
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7. 
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.90 
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.00 
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2 


.68 
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0. 
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.54 


6. 


.61 





.53 


4 


11 


0. 
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1470.11601 


-11 


.17 


8. 


.91 


41. 


,11 


5. 


.79 


0. 
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.751 


1545.08385 


27 


.02 


8. 


10 


16. 


.59 


4 


67 


0. 


179 





.906 


1545.86194 


-53 


.41 


8. 


.31 


1 


.59 


5. 


.47 


0. 


178 





.125 


1611.99999 


24 


.36 


16. 


.14 


-10 


.80 


22. 


01 


0. 


160 





.788 


1698.75509 


-64. 


.08 


9. 


61 


-22 


.52 


12. 


11 


0. 
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.268 


1699.76307 


16. 


.59 


8. 


.07 


-10 


.35 


7. 


.63 


0. 
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.553 



Note. — Note that for the iodine-free template exposures we do not measure 
the RV but do measure the BS and S index. Such template exposures can be 
distinguished by the missing RV value. 

a Barycentric Julian Date calculated directly from UTC, without correction 
for leap seconds. 

The zero-point of these velocities is arbitrary. An overall offset 7 rB i fitted to 
these velocities in Section |3.3| has not been subtracted. 

c Internal errors excluding the component of astrophysical jitter considered in 
Scction |3.3l 

Chromospheric activity index. 
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TABLE 4 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-40. 



BJD a RV b o- RV c BS o-bs S d Phase 
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0. 
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21 
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0. 
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.102 


1467.92695 


17. 


59 


2. 


.93 


-9. 


07 


6. 


.74 


0. 


147 





.542 


1468.99826 


57. 


77 


3. 


.09 


2. 


18 


2. 


.51 


0. 


145 





.782 


1469.73453 


28. 


31 


3. 


.45 


-11. 


76 


3. 


.57 


0. 


141 





.948 


1486.91164 


60. 


16 


3. 


.33 


-2. 


67 


6. 


.82 


0. 


145 





.801 


1521.81745 


35. 


22 


3. 


.96 


23. 


93 


6. 


.53 


0. 


143 





.633 


1698.11353 


-49. 


91 


3. 


.80 


4. 


39 


3. 


.60 


0. 


146 





.185 


1699.11678 


-31. 


14 


3. 


.84 


2. 


58 


4. 


.61 


0. 


150 





.410 


1701.07726 


16. 


71 


3. 


15 


-10. 


44 


1. 


.24 


0. 


145 





.850 


1853.75742 


-38. 


93 


2. 


73 










0. 


148 





.105 



Note. — Note that for the iodine-free template exposures we do not measure 
the RV but do measure the BS and S index. Such template exposures can be 
distinguished by the missing RV value. 

a Baryccntric Julian Date calculated directly from UTC, without correction 
for leap seconds. 

The zero-point of these velocities is arbitrary. An overall offset 7 ro i fitted 
to these velocities in Section |3.3l has not been subtracted. 

c Internal errors excluding the component of astrophysical jitter considered in 
Scction |3.3l 

Chromosphcric activity index. 



TABLE 5 

Relative radial velocities, bisector spans, and activity index measurements of 

HAT-P-41. 



BJD a 


RV b 


ctrv c 


BS 


CBS 




s d 


Phase 


Instrument 


(2,454,000+) 


(ms" 1 ) 


(ms -1 ) 


(ms - 1 ) 


(m s ) 












1375.97050 


92.63 


9.38 


25.18 


6.00 


0. 


.155 





.546 


Keck 


1375.97664 






20.47 


5.29 


0. 


155 


0. 


.549 


Keck 


1378.12394 


-78.22 


9.40 


0.84 


7.92 





.148 





.346 


Keck 


1378.58766 


52.31 


40.10 










0. 


.518 


FIESm 


1379.06917 


132.24 


9.90 


32.40 


7.80 





.153 





.697 


Keck 


1379.54578 


21.11 


51.30 













.874 


FIESm 


1380.47420 


-106.99 


49.00 













.218 


FIESm 


1381.10708 


-39.66 


10.23 


-1.39 


4.57 


0. 


.150 





.453 


Keck 


1381.61687 


112.71 


49.90 













.642 


FIESm 


1383.49357 


-79.79 


75.00 













.339 


FIESm 


1400.84669 


49.04 


12.32 


-60.12 


7.44 


0. 


139 





.780 


Keck 


1404.80595 


-52.53 


10.18 


-32.58 


7.48 


0. 


.148 





.250 


Keck 


1427.51377 


112.18 


40.10 










0. 


.679 


FIESh 


1428.54090 


-22.42 


91.50 













.060 


FIESh 


1429.38087 


-79.42 


50.00 













.372 


FIESh 


1431.54013 


-153.82 


68.60 













.173 


FIESh 


1465.87835 


5.42 


11.42 


-9.93 


6.17 





.145 





.919 


Keck 


1467.87264 


65.94 


10.34 


-11.95 


6.68 


0. 


.147 





.659 


Keck 


1469.89585 


-68.65 


10.18 


17.06 


4.94 


0. 


.143 





.410 


Keck 


1490.75265 


-113.04 


10.38 


9.02 


9.83 





.145 





.152 


Keck 


1500.73246 


24.76 


9.14 


12.08 


8.78 


0. 


.143 





.857 


Keck 


1704.09441 


-67.20 


10.27 


-1.09 


5.26 


0. 


.155 





.342 


Keck 


1814.94017 


16.90 


10.72 






0. 


.148 





.487 


Keck 



Note. — Note that for the iodine-free template exposures we do not measure the RV but 
do measure the BS and S index. Such template exposures can be distinguished by the missing 
RV value. 

a Barycentric Julian Date calculated directly from UTC, without correction for leap seconds. 

The zero-point of these velocities is arbitrary. An overall offset 7 rc i fitted to these velocities 
in Section l3.3l has not been subtracted. 

c Internal errors excluding the component of astrophysical jitter considered in Section |3.3l 

Chromosphcric activity index. 
e We indicate separately observations obtained with FIES using the medium-resolution fiber, 
and observations obtained with FIES using the high-resolution fiber. 
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Fig. 2. — Top panel: Kcck/HIRES RV measurements for 
HAT-P-39 shown as a function of orbital p has e, along with our 
best-fit circular model (solid line; see Table [Til l . and our best-fit 
eccentric model (dashed line). Zero phase corresponds to the time 
of mid-transit. The center-of-mass velocity has been subtracted. 
Second panel: Velocity O — C residuals from the best fit. The er- 
ror bars include a component from astrophysical jitter (43.0ms -1 ) 
added in quadrature to the formal errors (see Section l3.3l l. Third 
panel: Bisector spans (BS), with the mean value subtracted. The 
measurement from the template spectrum is included. The BS un- 
certainties are internal errors determined for each spectrum from 
the scatter of the individual BS values measured on separate or- 
ders of the spectrum; they do not include the unknown contribu- 
tion from stellar jitter. Fourth panel: Full width at half maxi- 
mum (FWHM) of the cross-correlation functions computed from 
the blue regions of the Keck/HIRES spectra, with the mean value 
subtracted. Bottom panel: Chromospheric activity index S. Note 
the different vertical scales of the panels. Observations shown twice 
are represented with open symbols. 

Based on these observations we measure the H- and 
L'-band magnitudes of the neighbor relative to HAT-P- 
41 to be AH = 2.46 ±0.06 mag and AL' = 2.6±2.0mag. 
Assuming the star is a physical companion to HAT- 
P-41, these magnitude differences are consistent with 
the neighbor being a ~ 0.7 M© K-dwarf, or roughly 
~ 3.5 mag fainter than HAT-P-41 in i band. 

The KeplerCam observations of HAT-P-41 described in 
Section I2T41 show HAT-P-41 to have an elongation in the 
wing of its PSF due to the companion, but the seeing is 
not good enough to resolve the stars given the magnitude 
difference. We use th e DAOPHOT and ALLSTAR programs 
(|Stetsonlfl987l 119921 ) to obtain PSF-fitting photometry 



Fig. 3.— Keck/HIRES observations of HAT-P-40. The panels 
are as in Figu re [21 The parameters used in the best-fit model are 
given in Tabic [TT1 



TABLE 6 

High-precision differential photometry of HAT-P-39. 



BJD a 


Mag b 


""Mag 


Mag(orig) c Filter 


(2,400,000+) 






54921.62065 


0.00694 


0.00127 


11.10550 i 


54921.62127 


0.00597 


0.00127 


11.10530 i 


54921.62187 


0.00631 


0.00126 


11.10640 i 


54921.62248 


0.00358 


0.00127 


11.10490 i 


54921.62311 


0.00600 


0.00127 


11.10610 i 


54921.62390 


0.00762 


0.00127 


11.10910 i 


54921.62452 


0.00787 


0.00126 


11.10670 i 


54921.62512 


0.00663 


0.00126 


11.10860 i 


54921.62574 


0.01035 


0.00126 


11.10940 i 


54921.62637 


0.01050 


0.00126 


11.11030 i 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a Barycentric Julian Date calculated directly from UTC, without 
correction for leap seconds. 

The out-of-transit level has been subtracted. These magnitudes 
have been subjected to the EPD and TFA procedures, carried out 
simultaneously with the transit fit. 

c Raw magnitude values without application of the EPD and TFA 
procedures. 
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Fig. 4.— Keck/HIRES and FIES/NOT observations of HAT-P- 
41. The panels are as in Figure[2] Here we use filled circles to show 
Keck/HIRES observations, filled triangles to show FIES/NOT ob- 
servations with the high resolution fiber, and open triangles to 
show FIES/NOT observations with the low resolution fiber. BS 
and S index measurements are only available for the Keck/HIRES 
observat ions . The parameters used in the best-fit model are given 
in Table [TT1 

for the two stars, and find that the neighbor is fainter in 
i band by Ai « 3.3 mag, with an uncertainty of at least 
0.1 mag. We conclude that the broad-band photometry is 
consistent with HAT-P-41 having a K-dwarf binary com- 
panion. At the distance of HAT-P-41, the 3'.'5 angular 
separation corresponds to a projected physical separation 
of ~ 1000 AU between the two stars. 

3. ANALYSIS 

3.1. Properties of the parent star 

We measured the stellar atmospheric parameters 
for each star using the Keck/HIRES Iodine- free tem- 
plate spectr a, together with the Sp ectroscopy Made 
Easy (SME; iValenti fc Piskunovl l l996| ) package, and the 
iValenti fc Fischer! (|2005l ) atomic line database. For each 
star, we obtained the following initial values and uncer- 
tainties: 

• HAT-P-39 - effective temperature T eSit = 6325 ± 
100K, metallicity [Fe/H] = 0.14 ± 0.1 dex, stellar 
surface gravity fogg* = 4.04 ± 0.1 (cgs), and pro- 




Hiif 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 
Time from transit center (days) 

Fig. 5. — Unbinned transit light curves for HAT-P-39, acquired 
with KeplerCam at the FLWO 1.2 m telescope. The li ght c urves 
have been EPD- and TFA-processed, as described in § 13.31 The 
dates of the events are indicated. Curves after the first are dis- 
placed vertically for clar ity. Our best fit from the global modeling 
described in Section l3.3l is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. The error bars represent the photon and background shot 
noise, plus the readout noise. 



TABLE 7 

High-precision differential photometry of HAT-P-40. 



BJD a 


Mag b 


"Mag 


Mag(orig) c Filter 


(2,400,000+) 




55456.60581 


0.00871 


0.00080 


10.22060 i 


55456.60739 


0.00713 


0.00067 


10.21900 i 


55456.60830 


0.00854 


0.00066 


10.22070 i 


55456.60930 


0.00851 


0.00066 


10.22100 i 


55456.61028 


0.00935 


0.00065 


10.22120 i 


55456.61124 


0.00789 


0.00065 


10.21980 t 


55456.61224 


0.00868 


0.00065 


10.22070 i 


55456.61319 


0.00695 


0.00065 


10.21870 i 


55456.61416 


0.00727 


0.00065 


10.21880 i 


55456.61513 


0.00996 


0.00064 


10.22200 t 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a Baryccntric Julian Date calculated directly from UTC, without 
correction for leap seconds. 

The out-of-transit level has been subtracted. These magnitudes 
have been subjected to the EPD and TFA procedures, carried out 
simultaneously with the transit fit. 

c Raw magnitude values without application of the EPD and TFA 
procedures. 
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2010 Sep 16 (FLWO) 



2010 Oct 17 (FTN) 



2011 Aug 16 (FTN) 



2011 Oct 13 (BOS) 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 
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Fig. 6. — Similar to Figure [5] here we show the follow-up light 
curves for HAT-P-40. The facility used for each each light curve is 
indicated next to the date of the event. 



TABLE 8 

High-precision differential photometry of HAT-P-41. 



BJD a 


Mag b 


°"Mag 


Mag(orig) c Filter 


(2,400,000+) 






55344.79616 


0.01161 


0.00077 


10.01060 i 


55344.79679 


0.00980 


0.00078 


10.00870 i 


55344.79744 


0.00586 


0.00078 


10.00480 t 


55344.79878 


0.00847 


0.00078 


10.00830 i 


55344.79946 


0.01379 


0.00078 


10.01310 i 


55344.80015 


0.01168 


0.00078 


10.01140 t 


55344.80083 


0.00987 


0.00077 


10.00820 i 


55344.80151 


0.01021 


0.00077 


10.00830 i 


55344.80219 


0.01031 


0.00076 


10.00870 i 


55344.80288 


0.01025 


0.00077 


10.00800 i 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a Baryccntric Julian Date calculated directly from UTC, without 
correction for leap seconds. 

The out-of-transit level has been subtracted. These magnitudes 
have been subjected to the EPD and TFA procedures, carried out 
simultaneously with the transit fit. 

c Raw magnitude values without application of the EPD and TFA 
procedures. 



5 arcsec 




2010 May 27 (FLWO) 

^ 

••'«•■.. ' 201 

2011 Jun 1 9 (FLWO) 



201 1 May 31 (FLWO) 




2011 Jul24(B0S) 




-0.15 -0.1 -0.05 0.05 0.1 0.15 

Time from transit center (days) 

Fig. 7. — Similar to Figure [5] here we show the follow-up light 
curves for HAT-P-41. The facility used for each each light curve is 
indicated next to the date of the event. 




Fig. 8.— //-band AO image of HAT-P-41 obtained with 
MMT/Clio2 showing the ~ 3'/5 neighbor to the south. 

jected rotational velocity v sin i = 12.7±0.5kms _1 . 

• HAT-P-40 - effective temperature T eff * = 6140 ± 
100 K, metallicity [Fe/H] = 0.25 ± 0.1 dex, stellar 
surface gravity logg* = 4.04 ± 0.1 (cgs), and pro- 
jected rotational velocity wsini = 6.7± O.Skms . 

• HAT-P-41 - effective temperature T cSir = 6007 ± 
100 K, metallicity [Fe/H] = 0.06 ±0.1 dex, stellar 
surface gravity log g+ = 3.68 ± 0.06 (cgs), and pro- 
jected rotational velocity v sin i = 20.6±0.5kms~ 1 . 



As described in our previous papers (e.g. i Bakos et al.l 
2010), these initial values were used to determine the 
qua dratic l imb-d arkening coefficients for each star from 
the iClaretJ (|2004[ ) tables. We then used the mean stellar 
density, determined from the normalized semimajor axis 
a/Ri,, together with the effective temperature and metal- 
licity to determine an initial estimate of the mass and 
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radius of each s tar from the Yonsei-Yale (YY) isochrones 
(|Yi et al.ll2001l) . This provided a refined estimate of the 
stellar surface gravity, which we fixed in a second itera- 
tion of SME for each star. For each system a third itera- 
tion did not change log g+ appreciably, so we adopted the 
values from the second iteration as the final spectroscopic 
parameters for each star. These parameters are listed in 
Table [9l In this same table we also list the available 
broad-band photometric magnitudes from the literature, 
and physical parameters, such as the stellar masses and 
radii, which are determined from the spectroscopic pa- 
rameters together with the model isochrones. As dis- 
cussed in Section |3~31 we adopt the parameters assuming 
a circular orbit for each planet. Some of the parameters, 
especially the derived stellar masses and radii, depend 
on the eccentricity; Table [10] lists the values for these 
parameters when the eccentricity is allowed to vary. 

The inferred location of each star in a diagram of a/i?* 
versus Teg*, analogous to the classical H-R diagram, is 
shown in Figure [9] In each case the stellar properties 
and their la and 2a confidence ellipsoids are displayed 
against the backdrop of model isochrones for a range of 
ages, and the appropriate stellar metallicity. For com- 
parison, the locations implied by the initial SME results 
are also shown (in each case with a triangle). 

We determine the distance and extinction to each 
star by comparing the J, H and K$ magnit udes from 
the 2MASS Catalogue (jSkrutskie et al.l I2006T >. and the 
V and Ig magnitudes fro m the TASS Mark IV Cata- 
logue (jDroege et alJ 12006ft . to the expected magnitudes 
fro m the stel l ar mo dels. We use the transformations 
by iCarpenterl ()2001[) to convert the 2MASS magnitudes 
to the photometric sy s tem o f the models (ESO), and 
use the iCardelli et al.l (|1989t ) extinction law, assuming 
a total-to-selective extinction ratio of Ry = 3.1, to re- 
late the extinction in each band-pass to the y-band ex- 
tinction Ay. The resulting Ay and distance measure- 
ments are given in Table [HJ We find total V band ex- 
tinctions of A v = 0.171 ± 0.135, 0.353 ± 0.127, and 
0.248 ± 0.134 mag for HAT-P-39 through HAT-P-41, re- 
spectively. For comparison, the total line of sight extinc- 
tions i n each direction, estimated from the lSchlegel et all 
(1998) dust maps, are: 0.146 mag, 0.7 24 mag, and 
0.513mag. Following lBonifacio et all pOOO) , wc estimate 
the expected distance-corrected extinction to each source 
to be 0.113mag, 0.306 mag, and 0.157mag, respectively. 
For HAT-P-39 and HAT-P-40 the measured and ex- 
pected values are consistent. For HAT-P-41 the broad- 
band photometry appears to point to a slightly redder 
star than expected based on the spectroscopic tempera- 
ture and expected extinction. As noted in Section |2.5[ 
HAT-P-41 has a close companion which is unresolved 
in the 2MASS or TASS catalogs, this companion is the 
probable cause of the discrepancy between the expected 
and observed magnitudes of HAT-P-41. 

3.2. Excluding blend scenarios 

The analyses of our reconnaissance spectroscopic ob- 
servations discussed in Section [2.2l rule out the most ob- 
vious astrophysical false positive scenarios for HAT-P-39 
through HAT-P-41. Additionally the spectral-line bisec- 
tor span (BS) analyses which we conducted (Figures [5J- 
[4]) provide constraints on more subtle b l end s cenarios 
similar to that presented in iTorres et al.l (|2004D . How- 



HAT-P-39 



4.0 




T3 

8.0 - 




Effective temperature [K] 
HAT-P-40 




120 - 



6500 6000 5500 

Effective temperature [K] 
HAT-P-41 



3.0 



8.0 









\C < 


• -•-■=>" ;;> / 


A 









Effective temperature [K] 

Fig. 9. — Model isochrones from IYi et alj II2001I ') for the metal- 
licities of HAT-P-39 (top), HAT-P-40 (center) and HAT-P-41 (bot- 
tom). For HAT-P-39 and HAT-P-41 the isochrones are shown for 
ages of 0.2 Gyr, 0.6 Gyr, and 1.0 Gyr to 3.0 Gyr in steps of 0.5 Gyr 
(left to right), while for HAT-P-40 the isochrones are shown for 
ages of 0.2 Gyr, 0.6 Gyr and 1.0 Gyr to 4.5 Gyr in steps of 0.5 Gyr 
(left to right). The adopted values of T e £f^ and a/Ft* are shown 
together with their la and 2rr confidence ellipsoids. In each plot 
the initial values of T e ff* and a/Ri, from the first SME and light 
curve analyses are represented with a triangle. 

ever, because HAT-P-39 and HAT-P-41 have high RV 
jitter, and consequently high BS scatter (~ 80ms _1 and 
~ 30ms _1 , respectively), relative to their RV semiampli- 
tudes (63.6 ± 10.4ms" 1 and 92.5 ± 11.6ms" 1 , respec- 
tively), the BS test provides a less stringent constraint 
on possible blend scenarios than it does in a case such 
as HAT-P-40 (BS scatter of - 10 ms" 1 and RV semi- 
amplitude of 58.1 ± 2.9ms" 1 ). To provide additional 
support for the planetary interpretation of the observa- 
tions of each system we conduct detailed blend analyses 
of the light curves (including both the HATNet discovery 
light curves, and all available f ollow-up light curves) and 
absolute photometry following lHartman et al.l (|201lD . 
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TABLE 9 

Adopted stellar parameters for HAT-P-39-HAT-P-41 assuming circular orbits 



HAT-P-39 HAT-P-40 HAT-P-41 

Parameter Value Value Value Source 

Spectroscopic properties 

T eSi , (K) 6430 ± 100 6080 ± 100 6390 ± 100 SME a 

[Fe/H] 0.19 ±0.10 0.22 ±0.10 0.21 ±0.10 SME 

vsini (kms^ 1 ) . . . 12.7 ±0.5 6.9±0.5 19.6 ± 0.5 SME 

« mac (kms- 1 ) . . . . 5.04 4.50 4.97 SME 

IW (kmr 1 ) . . . . 0.85 0.85 0.85 SME 

7 RV (kms" 1 ) 28.42± 0.28 -25.0 ±0.1 31.68 ± 0.61 DS or TRES 

Photometric properties 

V (mag) 12.422 11.699 11.087 TASS 

V-Ic (mag) 0.58 ±0.15 0.77 ±0.12 0.63 ±0.10 TASS 

J (mag) 11.424 ± 0.020 10.367 ± 0.023 10.006 ± 0.027 2MASS 

//(mag) 11.184 ±0.022 10.085 ± 0.018 9.777 ± 0.032 2MASS 

K s (mag) 11.157 ± 0.020 10.009 ± 0.018 9.728 ± 0.029 2MASS 

Derived properties 

M* (Mq) 1.404 ±0.051 1.512±o;?o9 1.418 ± 0.047 YY+a//?*+SME b 

R* (Rq) 1.625l°'°62 2.206 ±0.061 1.683^°-°^ YY+a/^+SME 

logg* (cgs) 4.16 ±0.03 3.93 ± 0.02 4.14 ±0.02 YY+a/i?*+SME 

L t (L Q ) 4.05 ±0.48 6.00 ±0.61 4.25 ± 0.41 YY+a/fl 4 +SME 

M v (mag) 3.21±0.14 2.83±0.13 3.16±0.12 YY+a/B t +SME 

Mx (mag,ESO) . . 2.12 ±0.10 1.49 ± 0.07 2.04 ± 0.06 YY+a/i?,+SME 

Age (Gyr) 2.0 ± 0.4 2.7±°'f 2.2 ± 0.4 YY+a/iJ*+SME 

A v (mag) c 0.171 ± 0.135 0.353 ± 0.127 0.248 ± 0.134 YY+a/i?* +SME 

Distance (pc) 642 ± 29 501 ± 16 344±g 2 YY+a/i?*+SME 

logR^ K d -4.85 ±0.07 -5.12 ±0.16 -5.04 ±0.04 Keck/HIRES 

a SME — "Spectroscopy Made Easy" package for the analysis of high-resolution spectra 
(Valcnti & Piskunov 1996). These parameters rely primarily on SME, but have a small dependence 
also on the iterative analysis incorporating the isochrone search and global modeling of the data, as 
described in the text. 

b YY+a/i?„+SME = Based on the YY isochrones lYi et a.1.1120011 1 . a/R t as a luminosity indicator, 
and the SME results. 

c V band extinction determined by comparing the measured 2MASS and TASS photometry for the star 
to the expected magnitudes from the YY+a/i?*+SME model for the star. We use the lCardelli et al.l 
( 1989) extinction law. 

Chromospheric activity index defined in Novcs et al. ( 1984) determined from the Keck/HIRES spectra 
following Isaacson & Fischer (2010). In each case we give the average value and the standard deviation 
from the individual spectra. 



TABLE 10 

Derived stellar parameters for HAT-P-39-HAT-P-41 allowing eccentric orbits' 1 



HAT-P-39 HAT-P-40 HAT-P-41 

Parameter Value Value Value Source 

M* (M Q ) 1.400io;o6g 1.504 ± 0.103 1.405 ±0.066 YY+a/K*+SME 

R* (Rq) 1.622ig'^ 2.422 ±0.154 1.525±g YY+a/R 4 +SME 

logg* (cgs) 4.16 ±6.10 3.85 ± 0.04 4.22 ± 0.07 YY+a/R t +SME 

L* (Lq) 4,0llJ'^ 6. 98±Q 97 3.68±J'°° YY+a/^+SME 

M v (mag) 3.22 ±0.31 2.67 ±0.18 3.31 ± 0.25 YY+a/i?*+SME 

M K (mag,ESO) . . 2.12 ±0.29 1.31 ± 0.14 2.25 ± 0.22 YY+a/R*+SME 

Age (Gyr) 2-0±°'.l 2 - 9 -°o.l 1.5 ± 0.6 YY+a/^+SME 

A v (mag) Q.lb&toXri 0.317 ±0.126 0.332 ±0.134 YY+a/K t +SME 

Distance (pc) 641 ti36 5 548 ± 36 31lj;p YY+a/i?*+SME 



a Quantities and abbreviations arc as in Table [9] which gives our adopted values, determined 
assuming circular orbits. We do not list parameters that are independent of the eccentricity. 
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In Figurc[TU]we show, for each system, the histogram of 
Ax 2 values between the best-fit transiting planet model 
and the best-fit blend model for simulated data sets hav- 
ing the same noise properties as t he observed residuals 
from the best-fit blend model (see IHartman et al.l 120111 
for a more detailed discussion). In this same figure we 
also show the A% 2 difference between the best-fit models 
applied to the observations. We find that for HAT-P-39 
and HAT-P-41 we can reject blend scenarios involving 
combinations of three stars with greater than 3tr and 5a 
confidence, respectively, based solely on the photome- 
try. For HAT-P-39 the detailed shape of the transit, as 
determined from the follow-up light curves, contributes 
most of the x 2 difference between the models, while for 
HAT-P-41 it is the lack of out-of-transit variations, as de- 
termined from the HATNet light curve, that contributes 
most of the x 2 difference. For HAT-P-40 we are unable 
to rule out blend scenarios based solely on the photome- 
try, however in this case the lack of BS or FWHM varia- 
tions rules out such blends. To quantify this, we simulate 
the cross-correlation function (CCF) of blended systems 
which could plausibly fit the photometric data (config- 
urations which cannot be rejected with greater than 5a 
confidence) and find that in all cases either the RV or 
FWHM of the blended configuration varies by several 
kms , or the BS varies by greater than ~ 100ms , 
greatly exceeding the observational limits on any such 
variations. Similarly for HAT-P-39 we find that stel- 
lar blend configurations which cannot be rejected with 
greater than 5a confidence predict greater than 500 m s _1 
variations in the RV or BS of the Keck spectra, which are 
well above the observational constraints. 

As discussed in Section |2~5| HAT-P-41 has a close 
neighbor which we estimate to be ~ 3.5 mag fainter in 
i band. While such a neighbor could, in principle, be 
eclipsed by an object that would produce a ~ 1% dip 
in the blended light curve, as we have shown here, the 
detailed shape of the light curve cannot be produced us- 
ing physically possible combinations of stars (i.e. stars 
with parameters determined from stellar evolution mod- 
els). We note that the KeplerCam observations show 
no evidence for variations in the flux centroid that cor- 
relate with the photometry, providing further evidence 
that the observed variation is not due to a deep eclipse 
in the poorly resolved neighbor. 

While we can rule out, for each of the systems, blend 
scenarios involving only stellar-mass components, we 
cannot rule out scenarios involving binary star systems 
with one component hosting a transiting planet. Indeed, 
for HAT-P-41 we find that including a ~ 0.7 Af© star in 
the system provides a slightly better fit to the photomet- 
ric observations. In this case we actually know that there 
is a faint companion (though it is unresolved in our light 
curves or the available absolute broad-band photometry 
measurements) , so it is reassuring to find that the blend 
analysis of the photometric data also points to the ex- 
istence of this companion. For HAT-P-39 we can rule 
out binary companion stars with M < 1.24 M© while 
for HAT-P-40 and HAT-P-41 we cannot rule out binary 
companions of any mass up to that of the mass inferred 
for the brighter star in the system. While massive binary 
companions in general should be easier to detect in the 
spectrum, if the two stars have very similar average ve- 




FlG. 10. — The histogram of A\ 2 values between the best-fit 
transiting planet model and the best-fit blend model for simulated 
data sets having the same noise properti es as the observed re sid- 
uals from the best-fit blend model (see Hartman ct al. 2011 for 
a more detailed discussion). The histogram is shown separately 
for each system. In each panel the arrow shows the x 2 difference 
between the models fitted to the actual observations. For HAT-P- 
41 the best-fit blend scenario is rejected based on the photometry 
with greater than 5ct confidence. For HAT-P-39 and HAT-P-40, 
blend scenarios which cannot be rejected based solely on the pho- 
tometry are rejected based on spectroscopic information (limits on 
variations in RV, and in the BS and FWHM of the spectral line 
profiles). 

locities, the resulting variations in BS, FWHM, and RV 
measurements can be less than the constraints set by the 
observations. 

We conclude that each system presented here contains 
a transiting planet, however we cannot definitively claim 
that these are all single stars. There is no evidence that 
either HAT-P-39 or HAT-P-40 is a binary system, so 
we treat each of these as single stars in the analysis that 
follows. For HAT-P-41 there is a resolved neighbor which 
we account for in our analysis of the system. 

3.3. Global modeling of the data 

We modelled the HATNet photometry, the follow-up 
photometry, and the high-precisio n RV measu r ement s 
using the procedure describ ed by IBakos et all (f 2010). 
Following the discussion by lEastman et al.l (|2012|) . we 
made two important changes t o our analysi s proc edure 
compared to what was do ne in IBakos et al.l (I2010D . As 
noted in section 3.5.3 of Eastman et al.l (|2012D . there 
is a common mistake in the implementation of the 
Metropolis-Hastings (M-H) algorithm for conducting a 
Markov Chain Monte Carlo (MCMC) analysis whereby 
the Markov chain is not increased when a proposed tran- 
sition is rejected. We discovered that the implementation 
we have been using for the analysis of HATNet planets 
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has made this mistake, and have corrected it for the anal- 
ysis of the planets presented in this paper. We found that 
this bug tends to inflate the errorbars on the determined 
parameters by a factor of a few parts in a hundred. Errors 
given in previous discovery papers may thus be slightly 
overestimated. The second significant change that we 
have made is to use ^fecosui and -y/esincj as jump pa- 
rameters, rather than ecosw and esinw. Previously we 
had been assuming uniform priors on the latter jump 
parameters, which amounts to assuming a linear prior 
on the eccentricity e, creating a bias towards measuring 
nonzero eccentricities. As other authors have noted, us- 
ing -y/ecosu; and y/esmuj leads to a uniform prior on e. 
We found that this change had a much more significant 
impact on our determined parameters than did correcting 
the bug in our implementation of M-H. For example, for 
HAT-P-39b the eccentricity that we find is 0.094 ± 0.086 
compared with 0.161 ± 0.094 when using our old jump 
parameters. 

We also made a few minor changes to the analysis 
for the particular planets presented in this paper. For 
the analysis of HAT-P-41 we allowed independent RV 
zero-points for the Keck/HIRES RVs, the high- resolution 
NOT/FIES RVs, and the medium-resolution NOT/FIES 
RVs. To account for the contribution from the neighbor 
to the photometry, we also fixed the third light in the 
i band to 4%, based on our PSF-fitting analysis of the 
KeplerCam observations of this system. 

For each planet we performed the fit both allowing ec- 
centricity to vary and fixing it to zero. The resulting pa- 
rameters for each system are listed in Table [Til assuming 
circular o rbits, and in Tabl e [T2l al lowing eccentric orbits. 
We use a lLucv fc Sweeney] (|197lD test to determine the 
significance of the measured eccentricities of each sys- 
tem. We find that the observations of all three systems 
are consistent with the planets being on circular orbits 
(the circular orbit hypothesis is rejected with 45% confi- 
dence, 86% confidence, and 91% confidence, for HAT-P- 
39b, HAT-P-40b and HAT-P-41b, respectively; for refer- 
ence the confidence level would need to be greater than 
99.7% for the detection of a nonzero eccentricity to be 
significant at the 3tr level, which is generally taken as 
a minimum level of significance). In the past we have 
generally presented parameters for systems allowing the 
eccentricity to vary, even in cases where the observations 
are consistent with a circular orbit, on the grounds that 
this provides a more conservativ e estimate of the errors . 
However, as discussed recently by I Anderson et al](|2012[) . 
the best-fit parameters which result from allowing the ec- 
centricity to vary are often biased relative to the circular 
orbit values, and in most cases further follow-up observa- 
tions, such as occultation observations, reveal the plan- 
ets to be on circular orbits after all. We therefore sug- 
gest adopting the circular orbit parameters as the most 
probable values for each planet. These are given in Ta- 
ble [TTJ For reference Table IT2T, lists the resulting param- 
eters when the eccentricities are varied. 

4. DISCUSSION 

We have presented the discovery of three new tran- 
siting planets which we show on mass-radius and equi- 
librium temperature-radius diagrams in Figure 1111 As 
seen in the mass-radius diagram planets generally have 
radii with 0.6 i?j < R < 1.5 i?j over a broad mass-range 



spanning over two orders of magnitude, except for in the 
range 0.4 Mj < M < 1.5 Mj where planets are found 
with radii as large as — 2.0 Rj (if WASP-12b is ex- 
cluded, then the mass range is 0.4 Mj < M < l.OMj) 
Applying the Kolmogorov-Smirnov test (e.g. I Press et al.l 
1992), we find that there is only a 0.7% chance that the 
masses of planets with M > 0.4 Mj and R > 1.5 Rj are 
drawn from the same distribution as the masses of plan- 
ets with M > 0.4 Mj and R < 1.5 Rj. The three plan- 
ets presented here fall in the population of large radius 
(R > 1.5 -Rj), sub- Jupiter-mass planets which we refer to 
as highly inflated planets. 

As has been rep eatedly noted (the earliest reference 
being iGuillotl l2005h the radii of close-in gas-giant plan- 
ets are strongly correlated with the degree of irradiation 
(variously traced by the planet equilibrium temperature 
estimated by adopting a constant albedo, typically zero, 
for all planets, and making an assumption about the heat 
redistribution, or traced by the bolometric surface flux) . 
As is evident in Figure [TTJ the degree to which planets 
are inflated depends on their masses, with lower mass 
planets showing a stronger correlation between temper- 
atur e and radius. Thi s has also been previously noted 
(e.g. lEnoch e~ al. 2012) and has been taken as evidence 
for some theoretical models of the inflation process (e.g. 
iBatvgin et all I201U lLaughlin et~aH 1201 If ). The plan- 
ets presented here generally follow the established em- 
pirical trends, though they are somewhat more inflated 
than other planets with comparable equilibrium tem- 
peratures and semima j or ax es. The empirical relation 
given bv lEnoch etH] (|2012l ). which gives a prediction 
for the radius as a function of T eq and a for planets with 
0.5 Mj < M p < 2.0 Mj, yields radii of 1.52 Rj, 1.63 Rj, 
and 1.58 Rj for HAT-P-39b th rough HAT-P-41b respec- 
tively. The formula given bv IBekv et all (|2011[ ). which 
uses T eq and [Fe/H] as indepedent variables and was 
derived for planets with 0.3 Mj < M p < 0.8 Mj pre- 
dicts radii of 1.31 i?j, 1.30 Rj and 1.37 Rj. In all cases 
the predicted radii are smaller than the me asured values 
(T57i?j, 1.73 Rj, and 1.68 Rj), though the lEnoch et all 
(|2012ft relation, which was determined including more 
recent discoveries, gives values that are closer to the ob- 
servations. 

We find that HAT-P-41 has a neighbor which has 
near-IR photometry consistent with it being a 0.7 M 
star at the same distance as HAT-P-41. HAT-P-41 is 
thus one of number of hot Jupiter host sta rs with close 
visually-resolved neighbors (e.g. HP 189733. lBakos et all 
[20061 HAT-P-1, IBakos et alJ l2007t and many others), 
though it is not known what fraction of these are phys- 
ical companions. Such companions may be responsi- 
ble for drivi ng planets into close-in orbit s via the Kozai 
Mechanism (jFabrvckv fc Tremainel [2007). a hypothesis 
which has gained traction recently with the discovery 
that many close-in planets are on high obl iquity orbits 
(e.g. lTriaud eTaT]|2010l; lAlbrecht et afl[2012]) . but a com- 
plete survey to determine the frequency of companions 
to hot Jupiter host stars in a statistically robust way has 
not yet been published. 

Finally we note that all three of these planets are good 
tar gets for measuring the Rossiter-M cLaughlin effect (R- 
M: lRossiterlll924 : lMcLaughlinlll924D as they orbit bright 
stars with relatively high projected rotation velocities. 
For HAT-P-39b and HAT-P-41b the expected R-M semi- 
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TABLE 11 

Adopted orbital and planetary parameters for HAT-P-39b-HAT-P-4Tb assuming circular orbits 





HAT-P-39b 


HAT-P-40b 


HAT-P-41b 


Parameter 


Value 


Value 


Value 



Light curve parameters 

P (days) 3.543870 ± 0.000005 

(BJD) a 2455208.75049 ± 0.00041 



Ti4 (days) 

T12 = T 34 (days) a 

a/R* 

</R* b 

RpjR* 

b 2 

b = a cos i I Ri, 

i (deg) 

Limb-darkening coefficients 

ci,i (linear term) 

C2 , i (quadratic term) . . . 



0.1745 ±0.0017 
0.0178 ±0.0017 
6.74 ± 0.25 
12.75 ±0.07 
0.0993 ± 0.0025 

n 199 +0.069 
u - lzz -0.061 
34Q+0 085 
U.O*3_ 120 

87.0 ± 1.0 



ci, 

C2, 

ci. 

C2- 



ci,R 

C2,R 

RV parameters 

K (ms" 1 ) 

e 

RV jitter (ms" 1 )" 1 . 

Planetary parameters 
M p (Mj) 

/.'. (R.i) 

C(M P ,R P ) e 

Pp (gcm~ 3 ) 

logg P (cgs) 

a (AU) 

(K) f 



0.1787 
0.3812 



0.2335 
0.3871 



4.457243 ± 0.000010 
2455813.17584 ± 0.00054 
0.2557 ±0.0014 
0.0196 ± 0.0009 

5 on+0.06 

8.48 ± 0.04 
0.0807 ±0.0014 
030+ 045 
q 174 +0.091 
u - 1 '^-0.077 

88.3 ±0.9 



0.2190 
0.3650 



2.694047 ± 0.000004 
2454983.86167 ± 0.00107 
0.1704 ±0.0012 
0.0166 ± 0.0009 
5 44+0-09 

13.01 ± 0.06 
0.1028 ± 0.0016 
u.u^a_ 029 

n 999 +0.088 
u '^^-0.093 

87.7 ± 1.0 



0.1908 
0.3746 
0.2658 
0.3814 



63.6 ± 10.4 
(fixed) 
43.0 



0s 

(F) (10 9 ergs" 



L cm 2 ) 



0.599 ± 0.099 

1 t-yi +0.108 
± -°' 1 -0.081 

0.09 
0.19 ±0.04 
2.77 ±0.09 
0.0509 ± 0.0006 
1752 ± 43 
0.027 ±0.005 
2.13 ±0.21 



0.1987 
0.3646 



58.1 ± 2.9 
(fixed) 
6.2 



0.615 ±0.038 
1.730 ±0.062 
0.36 
0.15 ±0.01 
2.71 ± 0.03 
0.06081°;°™ 

1770 ± 33 
0.029 ± 0.002 
2.22 ±0.17 



92.5 ± 11.6 
(fixed) 
33.4 



0.800 ±0.102 
1.685_ 051 

0.10 
0.20 ±0.03 
2.84 ± 0.06 
0.0426 ± 0.0005 
1941 ± 38 
0.028 ± 0.004 
3.20 ±0.25 



a Reported times are in Barycentric Julian Date calculated directly from UTC, without correction for leap seconds. 
T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. T14: total transit duration, 
time between first to last contact; T12 — T34: ingress/egress time, time between first and second, or third and fourth 
contact. 

k Reciprocal of the half duration of the transit used as a jump parameter in our MCMC analysis in place of a/R*. It 
is related to a/R t by the expression C/-R* = a/-R*(27r(l + e sin u j))/(PVl - b 2 y/l - e 2 ) JBakos et alll2010l 1. 
c Values for a quadratic law, adopted from the tabulations by Claret (2004) according to the spectroscopic (SME) 
parameters listed in Table [9l 

^ Error term, cither astrophysical or instrumental in origin, added in quadrature to the formal Kcck/HIRES RV errors 
such that x 2 P cr degree of freedom is unity. For HAT-P-41 we did not add a jitter term to the FIES/NOT RV errors 
because the formal errors for these observations exceeded the scatter in the RV residuals. 
c Correlation coefficient between the planetary mass M p and radius R p . 

^ Planet equilibrium temperature averaged over the orbit, calculated assuming a Bond albedo of zero, and that flux is 
rcradiatcd from the full planet surface. 

g The Safronov number is given by & = \ (V csc /V orb ) 2 = {a/ R P )(M P /M^) (see I Hansen fc Barman 200711 . 



Incoming flux pcr unit surface area, averaged over the orbit. 

amplitude is over 100 ms -1 . For HAT-P-40b the signal 
amplitude is lower, but this is compensated by the long 
duration of the transit. The three stars are also posi- 
tioned closely below (HAT-P-40) and above (HAT-P-39 
and HAT-P-4 1) the 6250 K effec tive temperature thresh- 
old found by iWinn et al.l ((2010) to separate planets on 
orbits that are well aligned with the spin axes of their 
host star s from planets that a re on high obliquity orbits 
(see also lAlbrecht et al.ll2012f ) . 
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TABLE 12 

Orbital and planetary parameters for HAT-P-39b-HAT-P-41b allowing eccentric 

ORBITS a 



HAT-P-39b HAT-P-40b HAT-P-41b 

Parameter Value Value Value 



Light curve parameters 

a/R* 6.74 ±0.77 5.39 ± 0.28 5.97 ±0.52 

C/R* 12.76 ±0.08 8.49 ±0.04 12.99 ± 0.06 

i (deg) S7.1±!j;;g 88. 0±?;| 87.9 ± 0.9 

RV parameters 

if (ms- 1 ) 63.8 ± 10.6 55.4 ± 2.5 95.6 ± 10.4 

Vecosw 0.094 ± 0.189 0.041 ± 0.064 -0.189 ±0.100 

v/esino; -0.008 ±0.263 0.300lo?3* -°- 293 ioi20 

ecosu 0.022 ±0.074 0.012 ± 0.018 -0.068 ± 0.039 

esinw -0.001 ± 0.114 0.092 ± 0.051 -0.105±0.086 

e 0.094 ±0.086 0.095 ± 0.048 0.139 ± 0.063 

lo (deg) 188 ± 115 83 ± 42 237 ± 38 

RV jitter (ms^ 1 ) 42.7 4.1 27.5 

Secondary eclipse parameters 

T a (BJD) 2455189.310 ±0.168 2455819.895 ± 0.052 2454985.092 ± 0.068 

T s 14 0.1741 ± 0.0378 0.3048 ± 0.0300 0.1397 ± 0.0241 

T 3 ] 12 0.0176 ± 0.0084 0.0239 ± 0.0029 0.0135 ± 0.0027 

Planetary parameters 

M p (Afj) 0.596 ±0.099 0.584 ± 0.035 0.812 ± 0.094 

R p (Rj) 1.565±g"f|| 1.900 ± 0.127 1.529tg ^ 

C(Mp,R p ) 0.08' 0.03 0.32' 

p p (gcm- 3 ) 0.19±°-J° 0.11 ± 0.02 0.28±°'°? 

log 9p (cgs) 2.77±6.14 2.60± 0.06 2.93 ± 6.09 

a (AU) 0.0509l°;ooJg 0.0607 ± 0.0014 0.0424 ± 0.0007 

T oq (K) 175llg7 7 1843 ± 60 1879 ± 88 

6 0.027 ±0.006 0.025 ± 0.002 0.032 ± 0.004 

(F) (10 9 crgs- 1 cm- 2 ) 2.12±°;|| 2.60± 0.34 2.82+° H 

a Quantities and definitions are as in Table fTTl which gives our adopted values, determined assuming 
circular orbits. Here we do not list parameters that arc effectively independent of the eccentricity. 




FlG. 11. — (Left): Mass-radius diagram of TEPs. HAT-P-39b through HAT-P-41b are indicated. Jupiter and Saturn are marked with 
filled squares. (Right): Equilibrium temperature versus radius, the mass of each planet is indicated by the color of its symbol. 
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N154Hr), and the NOAO Gemini/Keck time-exchange 
program (G329Hr). This paper presents observations 
made with the Nordic Optical Telescope, operated on the 
island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias. This paper uses observations obtained with 
facilities of the Las Cumbres Observatory Global Tele- 
scope. The Byrne Observatory at Sedgwick (BOS) is op- 
erated by the Las Cumbres Observatory Global Telescope 
Network and is located at the Sedgwick Reserve, a part 
of the University of California Natural Reserve System. 
Data presented in this paper are based on observations 
obtained at the HAT station at the Submillimeter Ar- 
ray of SAO, and the HAT station at the Fred Lawrence 
Whipple Observatory of SAO. The authors wish to rec- 
ognize and acknowledge the very significant cultural role 
and reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are 
most fortunate to have the opportunity to conduct ob- 
servations from this mountain. 
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